Probing the infrared quark mass from highly excited baryons 
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We argue that three-quark excited states naturally group into quartets, split into two parity 
doublets, and that the mass splittings between these parity partners decrease higher up in the 
baryon spectrum. This decreasing mass difference can be used to probe the running quark mass in 
the mid-infrared power-law regime. A measurement of masses of high-partial wave A* resonances 
should be sufficient to unambiguously establish the approximate degeneracy. We test this concept 
with the first computation of excited high-j baryon masses in a chirally invariant quark model. 
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Quantum Chromodynamics (QCD) has been thor- 
oughly tested in high-energy physics, through hadron 
jets, Drcll-Yan processes, e~e + annihilation, and deep 
inelastic scattering. Low-energy QCD manifests chiral 
symmetry breaking (%SB) that enhances quark masses in 
the infrared (IR), generating most of the visible mass in 
our universe, and that removes the degeneracy from the 
ground states of light hadron spectra (so light hadrons 
do not come in parity doublets). The excited light- 
quark baryon spectrum, with masses above that of the 
A-resonance (1232 MeV), has not been accessible to any 
form of perturbation theory, effective or fundamental. 

Inscnsitivity to chiral symmetry breaking, recently 
stressed by Glozman [H, 0, y] , and retrospectively present 
in the excited meson spectra computed in chirally invari- 
ant quark models 0, [J has led many hadron physicists to 
accept that spontaneous ;ySB, the salient feature of the 
low hadron spectrum, actually becomes less important 
for highly excited resonances. For these, chiral symme- 
try is asymptotically realized in Wigner mode, arranging 
hadrons in degenerate chiral multiplets @, 0] • We have 
realized that a formalization of this statement is that 
the ratio of the quark mass to its momentum (^) pro- 
vides a new perturbative parameter to study some as- 
pects of the excited hadron spectrum when the typical 
momentum becomes higher than the running mass. Rea- 
soning in reverse, we propose employing a cancellation 
in mass differences of very excited resonances, to access 
this small parameter, and hence the quark mass. We 
also show analytically and numerically that three-quark 
states naturally group into quartets with two states of 
each parity. Diagonalizing the chiral charge in terms of 
quarks, the quartet is split into two parity doublets, and 
all mass splittings tend to decrease when going higher in 
the spectrum. For simplicity and to minimize the impact 
of molecular meson- nucleon configurations [f| , we study 
the maximum-spin excitations A* of the Delta baryon. 

The u or d quark mass Q (understood as taken from 
the full two-point function of the field theory), is theo- 
rized to vary from to ~ 1 — 6 MeV at a high, perturbative 
momentum k scale (typically 3> Aqcd ~ 210 MeV), to a 
constituent mass of circa 300 MeV at vanishing momen- 
tum k ~ AgcD- This IR enhancement of two orders of 
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FIG. 1: IR enhancement of the light quark mass due to spon- 
taneous xSB. Shown are quark masses in the main approaches 
to QCD, all multiplied by an arbitrary factor to match them 
at quark momentum norm k — > 0. (References in the text.) 

magnitude is shown in Fig. [T] The simplest constituent 
quark model considers a constant mass [Io| . without 
the characteristic step-function of the Nambu and Jona- 
Lasinio model (also shown). Breaking chiral symmetry 
with a linear potential forces a power-law to — > -p- for 
m(k) <C k. Lattice QCD [ll| and Schwinger-Dyson equa- 
tions both Euclidian approaches to QCD, also show 
power-law decays of the light quark mass. However lit- 
tle has been said about this experimentally, for k ~ 1.5 
GeV, where we already have m(k) <C k. 

In Hamiltonian dynamics, the quark mass appears in 
the Dirac spinors U K \ that govern the quark-quark inter- 
action. In the heavy quark limit, spin-tensor potentials 
have been successfully derived by expanding in orders of 
k/m{k). For light quarks the opposite, ultrarelativistic 
limit of k ^> m(k) is relevant. Spinors are then conve- 
niently expanded in the inverse ratio m(k)/k, or 
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with E(k) = y/ k 2 + m(k) 2 , having kept the leading chi- 
rally invariant term, and a next order chiral symmetry 
breaking term. Non-chiral, spin-dependent poten- 
tials in the quark-quark interaction originate from the 
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second term in an expansion of ffQ CD in the weak 
sense (that is, not of the Hamiltonian operator itself, but 
a restriction thereof to the Hilbcrt space of highly excited 
resonances, where (k) is large): 



(n\H^ CD \n') ~ {n\H® CD \n') 

(2) 

To access the quark mass from experiment one can ex- 
ploit the approximate parity degeneracy, that would be 
exact in the presence of only the first term H® CD , of ex- 
cited baryons. This degeneracy follows from invariance 
under transformations generated by the chiral charge flij 
Ql = J dx^(x)~f 5 ^-ip(x) and [Q%,H] = 0. However 
Chiral Symmetry is spontaneously broken by the ground 
state, Q||0) ^ providing a large quark mass in the 
propagator, m(fc), pseudo-Goldstone bosons (tt, K, rf), 
and the loss of parity degeneracy in ground state baryons. 
Substituting the spinors, and in terms of Bogoliubov- 
rotated qq normal modes B and D, Q5 becomes 
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In the presence of Spontaneous x^B, m(fc) 7^ 0, and 
the last term realizes chiral symmetry non-linearly in 
the spectrum as it creates/destroys a pion. But when 
(k) m(k), the cr -k-term dominates and chiral symme- 
try is realized linearly (with only quark counting opera- 
tors flipping parity and spin). This happens for a highly 
excited baryon resonance whose constituents have a mo- 
mentum distribution peaked at higher momenta than the 
IR enhancement of m{k) (Fig. [TJ) . 

For the ground-state resonances in excited spin chan- 
nels, three-quark variational wavefunctions are relevant 
and a reasonable phenomenological guide. When Q5 acts 
on such wavefunction \N) — FijkB\B^B\\^l) , the re- 
sult also contains three quarks, but one of them is spin- 
rotated from Bk\ to a- k\\>Bk\>- Successive application 
of the chiral charge spin-rotates further quarks, changing 
each time the parity of the total wavefunction. However 
the sequence of states is closed since a-ka-k = I. In fact, 
starting with an arbitrary spin-symmetric wavefunction 
with parity P, one generates our claimed Baryon Quartet 
of equal isospin (isospin index dropped) : 
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which is a natural basis to discuss chiral symmetry 
restoration in baryons, through wavefunctions that are 
linear combinations \N) = X] c a \N a ). This reducible rep- 
resentation of Qg =3 corresponds to [(3/2,0) ® (0,3/2)] 2 
of SU(2)l x 5J7(2)ks and can be embedded in larger ones 
including nucleons [7j . 

Because [Qs,H] = 0, the two can be diagonalized si- 
multaneously The representation of the chiral charge in 
the (non-orthonormal) quartet coordinates breaks in two 
blocks since Q5 changes parity, or using the square charge 



Qt 



3 2 

6 7 



(5) 



Immediately one sees that the two linear combinations 
N0—N2 and No+3N 2 diagonalize the square chiral charge 
in the positive parity sector (if the original wavefunction 
had positive parity) , with Ni — N3 and 3iVi + N3 doing 
so for negative parity. The quartet then separates into 
two doublets connected by the chiral charge 

Q 5 (N -N 2 )=N l -N 3 Q 5 (N +3N 2 ) =3(3JVi+JV a ) 
Qs(iVi - N 3 ) = N - N 2 Q 6 (3JVi + JV 3 ) =Z{N Q + 3N{$>) 

Two doublets appear with different eigenvalues of Q\, 1 
and 9 respectively, and the interdoublet splitting becomes 
a dynamical question (we will argue shortly that it is 
small for highly excited baryons). However, the splitting 
within the doublet must vanish asymptotically. Even for 
fixed (not running) quark mass, when the typical kinetic 
energy is high enough (k) » m, the effects of the quark 
mass are negligible. Parity doubling then boils down to 
whether the interaction terms are also chiral symmetry 
violating or not. 

We now exploit the smallness of the \M P=+ — M p= ~\ 
mass difference with increasing angular momentum to 
guide experiment in obtaining information about the 
running quark mass. For this we establish how the j- 
dependence of the splitting is related to the fc-dependence 
of m(k). The M ± in our proposed study are the masses 
of the ground state quartets of the A spectrum, with par- 
ity ± and in the limit of large total angular momentum 
j >> 3/2. The two approximately degenerate masses 
M + and M~ are in the same leading linear Regge tra- 
jectory, phenomcnologically fixing their j-dependence to 
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The parity of the ground state alternates between + and 
— as the angular momentum steps up by one. Large j 
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FIG. 2: Typical momentum distributions of increasingly 
excited A 3 / 2 , . ..A 13 / 2 resonances overlap less and less 
with the dynamically generated IR quark mass so that 
\ri) becomes small. (Illustrative variational 
wavefunctions for a linear potential with string tension a = 
0.135 GeV 2 , not normalized for visibility). 

is equivalent to large quark orbital angular momentum 
/ 3> 1 since the spin is finite, and also to a large average 
linear momentum (k) 3> Aqcd- Figure [2] shows how the 
wavefunctions of high-j states overlap with the running 
m(k). In these states (m(k)/k) is a useful perturbative 
parameter. From the relativistic virial theorem (lBj be- 
tween kinetic and total energies, 



(k) -> c 2 M= 



C-2 

Ja 



'J 



(8) 



(where c 2 is a constant; for instance in the case of a linear 
potential model, and 3 quarks in a baryon, c-2 = 1/6). 

The first term in Eq. © cancels out in the difference 
\M+ - M~\ « M ± (while both M ± are dominated by 
the chirally invariant term, their difference stems from 
the dynamically generated quark mass) thus exposing the 
second term in cq.@, proportional to ( m [ ) , viz. 
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(the factorization is allowed by the mean-value theorem 
at the price of a constant that we do not attempt to deter- 
mine here). This equation is analogous to the renowned 
Gell-Mann-Oakes-Renner relation A/ 2 = —m q ^jP- but 
active when chiral symmetry is realized linearly, as in 
the high-baryon excitations we examine. 

Notice that H® CD ' differs from H® CD due to the - 
sign of the (er- k)j term in the (m(k)/k) expansion of the 
quark i spinor in eq. ((2|). This changes the sign of spin- 
dependent terms in H^ CD , led in high j by the spin- 
orbit term. In H® CD the spin-orbit is crucial to correct 
the angular momentum in the centrifugal barrier term 
from to the chirally invariant L^ 2 + 2L^ • Sj = J; 2 — |. 
But the opposite sign of the spin-orbit term of H^ CD 
produces a \M + — M~\ mass difference, enhanced in the 
IR by (m(k)/k). Since the centrifugal barrier scales in 



FIG. 3: Parity doubling in the spin-excited A spectrum for 
a string tension a = 0.135 GeV 2 and the mass gap angle 
matched to the lattice data. A qqq variational Monte-Carlo 
computation of the masses with a chiral Hamiltonian shows 
that the ground-states with parity ± for each j = 3/2 ... 13/2 
quickly degenerate. Experimentally the degeneracy can be 
claimed for the 9/2 states only. 

the large j limit like the baryon mass M ± , and the spin- 
orbit term comes with one power of j less than it, 
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Combining the j-power-law analysis of eqs. (O to (fTU]) . 
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that links the IR enhancement of the quark mass to 
baryon spectroscopy in a usable way. An experimental 
extraction proceeds by just fitting the exponent —i of j 
in the splitting \M + — M~\ cx j~ l . Then, we obtain, 



(-^L x y/j) oc and m{k) oc fc~ 2,;+2 . (12) 



The same exponent i in eq. (| 12[) can be obtained from 
the fit [I3| to the |M+ - M~| with increasing j\ Though 
the quark mass itself is gauge-dependent, our analysis 
suggests that its power-law exponent is not, since it is 
directly related to an observable. 

We next contribute the first competitive chirally in- 
variant quark model computation of parity doubling in 
the excited baryon spectrum, shown in figure [3) employ- 
ing a well-known model of Coulomb-QCD that can be 
seen as a field theory upgrade of the Cornell model, 

H = -g. J dx* t (x)a-A(x)*(x) +Tr J dx(E-E + B-B) 

+ J dx¥ q (-ia-V + f3m)* q ~[dxdypZV L (\x-y\) P ; 

with the kernel Vl{t) = or with string tension a = 0.135 
GeV 2 , coupled to the color charge density p a (x) = 
*t( x )T a *(x) + f abc A h (x) ■ U c (x) . We solve the BCS 
gap equation and employ both the so calculated m(k) 
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FIG. 4: All model quartet splittings decrease for high j. Left: 
intradoublet splitting for the first parity doublet. Right: in- 
terdoublet splittings between natural and unnatural parity 
states. The latter is smaller for low j. 



quark mass as well as lattice computation in Landau 
and Coulomb gauges [llj. The model has the same 
chiral structure of QCD[4j, satisfying the Gell-Mann- 
Oakes-Rcnncr relation and other low-energy theorems, 
and allowing computations of static Nir observables 17 1 . 
Pauli antisymmetrization of wavefunctions, spin summa- 
tions and the 9-dimensional integrations are performed 
in a computer program, employing Monte-Carlo integra- 
tion [181] . to compute the excited A* spectrum. We pro- 
ceed variationally and employ several types of wavefunc- 
tions. While a suboptimal variational basis for low-lying 
states, the quartet states coincide with ii-eigenstates 
high in the spectrum. We use both the quartet and the 
fixed L variational basis. 

As can be seen from figure 21 within our variational 
and Monte-Carlo integration errors, the model splittings 
drop with j, as predicted analytically; a model statement 
equivalent to the generic cq. j9]) reads 



M+ - M- = 3 
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If precise data becomes available at ELSA or Jefferson 
Lab for the Aj with J = 7/2, 9/2, 11/2, 13/2 •• • par- 
ity quartets, one should be able to distinguish between 
the 1/j fall of \M + — M~\ for non-chiral models (with 
a constant difference {M + — AI~) 2 between the Rcggc 
Trajectories [HJ), and the faster drops for chiral theories 
such as QCD. Since the two doublets are closely degener- 
ate, both positive and negative parity ground states will 
have a nearby resonance with identical quantum num- 
bers. Given their width, it is likely they will only be 
distinguished by very careful exclusive decay analysis. 

Further, if the A spectrum could be measured up to 
high resonance masses, beyond the IR m(k) enhance- 
ment, and a lattice calculation of (H® CD ) became avail- 
able, an almost direct measurement of the current quark 
mass would follow in a regime of logarithmic, rather than 
power-law running. We advocate that a measurement 



of masses of high partial-wave A resonances with an 
accuracy of 50 MeV should suffice to establish the ap- 
proximate degeneracy and test the concept of running 
quark mass in the IR. The presented concepts should 
also further motivate analysis of empirical datafljj such 
as EBAC (Excited Baryon Analysis Center) at Jefferson 
Lab as well as Lattice QCD collaborations who are ad- 
dressing the baryon spectrum. 
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